PROCEEDINGS 

OF SCIENCE 



NLO event samples for the LHC 



O M.V. GARZELLI, Adam KARDOS and Zoltan TROCSANYI*^ 

Institute of Physics, University of Debrecen, H-4010 Debrecen P.O.Box 105, Hungary 
^ E-mail: Z . TrocsanyiSatomki . hu 

o 

We introduce a twiki page with collections of generated Monte Carlo event samples in proton- 
proton collisions at LHC energies including a heavy quark-antiquark pair in the final state. These 
(— I samples are generated with the POWHEG method and can be used to prepare distributions at the 

NLO accuracy with first radiation treated according to the parton shower approach. Information 
I related to each event is stored in the form prescribed by the Les Houches Accords. Standard parton 

^ shower Monte Carlo programs can be used to further evolve these events, and simulate events at 



X 



the hadron level, ready for almost arbitrary experimental analysis. Currently the available final 
states are the following: (i) t + t, (ii) t + t + H, (iii) t + t + jet, while the generation of several other 
final states is in progress. 
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In recent years a lot of NLO QCD calculations have been presented in the literature. High- 
energy experiments, notably at the LHC, have and will be benefited by the progress in our com- 
putational ability to deal with higher order corrections in scattering amplitudes with many partons 
involved. In order though to get the optimum benefit and to produce predictions that can be di- 
rectly compared to experimental data at the hadron level, a matching with parton showers (PS) and 
hadronization is ultimately inevitable. 

In order to construct a generic interface between parton showers and matrix element calcu- 
lations at the NLO accuracy, we have chosen to combine the HELAC-NLO [3, 4] and POWHEG 
[1,2] approach respectively. Our goal is to apply our method to the hadroproduction of a tt pair 
in association with Standard Model boson(s) or fermion(s) (denoted X). Due to their largest mass, 
the t-quarks play an outstanding role in the Standard Model and top physics is expected to yield a 
lot of new experimental information at the LHC (see contributions in this proceedings). The first 
applications of our project were the processes pp — ^ tt -i- jet [5] and pp — tt + H [6]. In this con- 
tribution we introduce a web-page that is planned to contain event samples in the form of the Les 
Houches accord [7] for such processes. These events are readily usable for producing distributions 
at the hadron level including NLO QCD accuracy for the hard scattering process, and we expect 
that they will be very useful for experimental analyses. 

We performed our calculations using the POWHEG-BOX [8], which requires the following in- 
gredients: 1) flavor structures of the Born and real radiation subprocesses, 2) Bom-level phase 
space, 3) squared matrix elements with all incoming momenta for the Bom and the real-emission 
processes are built using amplitudes obtained from HELAC-Oneloop [9] and HELAC-PHEGAS 
[10], respectively. The matrix elements in the physical channels are obtained by crossing, 4) color- 
con^elated squared matrix elements ai^e taken from HELAC-Dipoles [11], 5) we use the polai^- 
ization vectors to project the helicity-correlated amplitudes to Lorentz basis for writing the spin- 
con^elated squared matrix elements. With this input POWHEG-BOX can be used to generate events 
at the Born level plus first radiation and store those in Les Houches Event Files (LHEF). Then one 
can choose any parton shower (PS) Monte Carlo program for generating events with hadrons. 

In the POWHEG-BOX the first emission is the hardest one measured by transverse momentum 
which is can also be chosen the ordering variable in PYTHIA. If the ordering variable in the shower 
is different from the transverse momentum of the parton splitting, such as in HERWIG, then the 
hardest emission is not necessarily the first one. In such cases HERWIG discards shower evolutions 
(vetoed shower) with larger transverse momentum in all splitting occuring after the first emission. 
In addition, a truncated shower simulating wide-angle soft emission before the first emission is 
also needed in principle, but its effect was found small [12, 5, 6]. As there is no implementation of 
truncated shower in HERWIG using external LHE event files, the effect of the truncated showers is 
absent from our predictions. 

As default, we check the consistency between real-emission. Born, color-coiTclated and spin- 
correlated matrix elements in randomly chosen phase space regions by taking the soft and coUinear 
limits of the real-emission matrix elements in all possible kinematically degenerate channels. We 
compare the virtual contributions to the predictions of other programs for computing one-loop 
amplitudes, such as MAD LOOP [13], if available. In all cases we found agreement up to 5-6 digits. 

The processes we considered were also studied extensively at the NLO accuracy in the litera- 
ture, which enables us to make detailed checks of our predictions. For each process we compared 
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Figure 1: Distributions of transverse momemtum (upper panels) and rapidity (lower panels) of the jet in the 
process pp tt + jet (left) and of the Higgs particle in the process pp — > tt + H (right). The insets below 
each plot show the ratio of the predictions (PowHel/NLO). 



the cross sections at the NLO accuracy to already published predictions, if available, and found 
agreement within the statistical uncertainty of our integrations. Finally, we also compared differ- 
ential distributions based on events already including first radiation from the POWHEG-BOX with 
the corresponding NLO predictions. As examples, we show transverse momentum and rapidity 
distributions in Fig. 1 for the relevant processes. The difference between the NLO and PowHel 
(POWHEG+HELAC) predictions are expected to be beyond NLO accuracy [1]. 

Next we present selected predictions 1 1 +X production with parton shower and hadronization 
effects at LHC. We used the event files from our webpage. Together with each set of events, one 
can download a sample analyis file to produce the plots presented here. In each plot we include 
our selection cuts, which can be altered almost freely to perform different analyses on the same 
events. Also the set of physical parameters used during the generation of the events can be found 
together with the event files. We used the last version of the PS programs, PYTHIA 6.425 [14] and 
HERWIG 6.520 [15], as well as for jet reconstruction. Fast Jet 2.4.3 [16]. 

In Fig. 2 we present distributions of the scalar- sum of transverse momenta in the event, H±for 
the two processes pp ^ tt + jet and pp ^ tt + H . The selection cuts implemented for these 
analyses are specified in Refs.[5] and [6], respectively. We can observe a softening of the spectra. 
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Figure 2: //x -distributions of the tt + jet (left) and of pp ^ tl + H (right) processes after decay of 

the heavy particles and after full shower Monte Carlo. The lower panels show the ratio of the predictions at 
the decay level to the shower MC level (PowHel+PYTHIA) and the ratio of the predictions at the hadron 
level using different shower MC codes (PowHel+HERWIG/PowHel+PYTHIA). 

when the shower is turned on, as expected on the basis of the parton shower approach. 

Interfacing NLO calculations, as structured in HELAC-NLO, with PS and hadronization ef- 
fects, within the POWHEG framework, open the exciting possibility of realistic, precise and reli- 
able simulations. We forsee large experimental and phenomenological potential of our approach. 
The event files produced by the POWHEG-BOX for processes pp — tt-i-X, where X is a hard ob- 
ject (SM boson, jet etc), together with the corresponding version of the program are available at 
http://grid.kfki.hu/twiki/bin/view/DbTheory/TthProd. 
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